Nonequilibrium growth of patchy-colloid networks on substrates 
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Patchy colloids with highly directional interactions are ideal building blocks to control the local 
arrangements resulting from their spontaneous self-organization. Here we propose their use, com- 
bined with substrates and nonequilibrium conditions, to obtain novel structures, different from those 
of equilibrium thermodynamic phases. Specifically, we investigate numerically the irreversible adhe- 
sion of three-patch colloids near attractive substrates, and analyse the fractal network of connected 
particles that is formed. The network density profile exhibits three distinct regimes, with different 
structural and scaling properties, which we characterize in detail. The adsorption of a mixture of 
three- and two-patch colloids is also considered. An optimal fraction of two-patch colloids is found 
where the total density of the film is maximized, in contrast to the equilibrium gel structures where 
a monotonic decrease of the density has been reported. 
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The past few years have witnessed a sustained interest 
in the self-organization of patchy colloids. These par- 
ticles, with functionalized surfaces, yield new features 
such as anisotropic interactions, control of the valence, 
and the formation of permanent electrical dipoles, paving 
the way to the development of novel materials with fine 
tuned mechanical, optical, and thermal properties [US]. 
Understanding how nonequilibrium conditions influence 
this self-organization is crucial to develop strategies to 
design new materials, as the novel structures emerge at 
very low temperatures, where thermal and mechanical 
equilibration might be difficult to achieve under normal 
experimental conditions. 

Theoretical and experimental studies of patchy colloids 
have been focused on their cooperative behavior in solu- 
tion [4- 7 , where several models were considered with 
the aim of describing a range of more complex build- 
ing blocks such as amphiphilic molecules, colloidal clays, 
proteins, and DNA nano-assemblies [TJ 0J |6j [SHE] • By 
contrast, the investigation of self-organization of the sim- 
plest of these models at planar substrates is only just 
beginning p~4j[15]. Theoretical studies have shown that, 
in the presence of a substrate, a rich equilibrium phase 
diagram emerges with very unusual properties such as, 
e.g., two wetting transitions and a nonmonotonic surface 
tension [15] . These works have been restricted to equilib- 
rium features; however, it is widely recognized that the 
use of substrates might improve the degree of control of 
aggregation, especially under nonequilibrium conditions 
P~6j[17]. For example, a growth direction can be denned 
which allows the use of kinetic features to control the film 
structure [T8] . 

In this Letter we address the nonequilibrium adsorp- 
tion of patchy colloids on substrates and characterize the 
resulting network of connected particles. In deep con- 
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FIG. 1. (color online) Density profile of the colloidal net- 
work on a substrate showing three different regimes: surface 
layer (green, left); liquid film (red, middle); interfacial re- 
gion (blue, right). Top: Snapshot of a typical configuration, 
where each stick corresponds to a colloid-colloid connection. 
Bottom: Density, p*(z), as a function of the height, z, after 
the adsorption of 40 particle layers, averaged over 500 inde- 
pendent realizations, pi is the density of the liquid film and 
Z1/2 is the film thickness, defined as the height at which the 
density is pi/2. 



trast with equilibrium films, a fractal network is assem- 
bled with a fractal dimension compatible with the one re- 
ported for Diffusion Limited Aggregation (DLA) [T91I2Q]. 
We systematically analyze the dependence of the network 
structure on the substrate size and diffusion coefficient of 
the colloids in solution (which can be controlled experi- 
mentally, for example, by the thermostat temperature). 
We show that, although the density of the film strongly 
depends on the diffusion coefficient, the fractal dimension 
is resilient over a wide range of growth conditions. 

It has been shown that the distribution of patches af- 
fects the aggregation process and, consequently, the equi- 
librium bulk structures [2j [5j [7] [I3j EEJ [22] . For exam- 
ple, control of the valence allows tuning up the density 
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FIG. 2. (color online) Right: Dependence of the density 
of the first adsorbed layer, p*(0), on D, for values ranging 
from 10 ~ 3 to 1, and substrates with lateral size, L, ranging 
from 16 to 128. Right: Snapshots of the first adsorbed layer 
(directly in contact with the substrate) for different diffusion 
coefficients, D, namely, 10 -3 , 10 -2 , 10 _1 , and 1, on a sub- 
strate with lateral size 32 in units of the particle diameter. 



and temperature of both the gas-liquid and sol-gel criti- 
cal points [5j[23][24]. Here, we consider three-patch col- 
loids and investigate the density profile of the adsorbed 
nonequilibrium network. We then proceed to investigate 
the adsorption of three- and two-patch colloids and char- 
acterize the dependence of the film density on the concen- 
tration of the mixture. Studies of the coexisting thermo- 
dynamic structures of these mixtures in solution reveal 
a monotonic decrease of the density with the fraction of 
two-patch colloids [22j |23j [25]. Somewhat surprisingly, 
we have found that the density of the adsorbed network 
film, under nonequilibrium conditions, increases with the 
concentration of two-patch particles before it decreases, 
exhibiting a well defined maximum at an intermediate 
concentration. 

In the following we quantify these results and illus- 
trate the full scaling behavior of the nonequilibrium ad- 
sorbed film. In order to simulate the colloidal adsorp- 
tion we propose and use a stochastic model based on the 
nonequilibrium Monte Carlo (MC) method. Patchy col- 
loids are frequently described as spherical particles with a 
short-range repulsive core and patch-patch attractive in- 
teractions. The attraction is truncated at a certain angle 
around the center of the patch [26]. Here, we introduce 
the concept of interaction range, as a region around the 
patch (on the surface of the colloid) where the interaction 
with other patches is effective. This region is typically 
much smaller than the colloid diameter. When a colloid 
approaches the patch of another one, the binding prob- 
ability corresponds to the probability of overlapping of 
their interaction ranges. Thus, our model contains a sin- 
gle control parameter (the interaction range) that can be 
tuned to mimic different patch-patch potentials. 

In the presence of a substrate, two characteristic 
timescales can be identified: one related to the flux of 
colloids towards the substrate (inter- arrival time) and the 
other to the binding between patches (binding time). At 
low concentration of colloids, the inter-arrival time can 
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FIG. 3. (color online) Data collapse for the liquid film and 
interfacial region, including the dependence on the substrate 
lateral size L and the diffusion coefficient D. Results ob- 
tained after adsorbing 60 particle layers and averaging over 
500 samples for L = {16,32}, 200 samples for L = 64, and 
100 samples for L = 128. 



be considered much larger than the binding time. There- 
fore, we assume that colloids arrive one at a time towards 
the substrate and adhere instantaneously. In addition, 
we consider chemical bonds between the patches, which 
are highly directional and assumed irreversible within the 
timescale of interest. In order to describe the colloid 
motion in solution, a Brownian algorithm is considered, 
adapting ideas previously implemented for Molecular Dy- 
namics (see, e.g. Ref. [5 ). In this limit, only three pro- 
cesses can occur, namely, collision with the surface, with 
a pre- adsorbed colloid, or with the solvent. At every col- 
lision with the solvent (Poisson process), a new velocity is 
randomly assigned, drawn from the Maxwell-Boltzmann 
distribution at the thermostat temperature. By defining 
the thermostat temperature and the properties of the 
Poisson process we can adjust the diffusion coefficient. 
When the collision is with a pre-adsorbed colloid and the 
contact point is within its interaction range, the possibil- 
ity of binding is checked based on the binding probabil- 
ity: if binding is successful, the position of the landing 
colloid is adjusted based on the patch-patch orientation; 
otherwise, an elastic collision occurs. The substrate is 
considered attractive and, therefore, collisions with the 
substrate always result in irreversible colloidal binding. 

Our model shares important features with Diffusion 
Limited Deposition (DLD), an extension of the famous 
DLA developed to account for the growth of films on 
fibers and surfaces [19]. In DLD, particles diffuse one 
after the other and adhere to the first pre-adsorbed par- 
ticle. For patchy colloids, since a link is established only 
when there is an effective overlap between interaction 
ranges, the colloid does not, necessarily, bind during the 
first collision. Besides, somewhat artificially, in the con- 
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FIG. 4. (color online) Snapshots of a region in the liquid 
film for different fractions of two-patch colloids, td> From 
left to right, top to bottom: 0, 0.3, 0.6, and 0.8. Snapshots 
obtained from the adsorption on a substrate with a lateral 
size of 32 particle diameters and 40 particle layers. Three- 
patch colloids are in blue (dark), two-patch colloids are in 
green (light), the (red) spheres on the surface of the colloids 
represent the patches, and the (red) sticks are the connections 
between colloids. 

tinuum version of DLD a fixed step length is considered 
(typically of the order of the diameter of the particle) 
and, consequently, the diffusion coefficient only affects 
the time scale. Here, we use a more realistic descrip- 
tion where the inter-collision time with the solvent, and 
consequently the step length, is a function of the thermo- 
stat temperature. As we will show below, an interesting 
dependence on the diffusion coefficient emerges. 

We start with an empty planar square substrate, with 
lateral size L, defined in units of the particle diameter, 
and assume periodic boundary conditions in the hori- 
zontal plane. Iteratively, spherical colloids, with three 
patches uniformly distributed on their surface, are re- 
leased one after the other and diffuse until they bind ei- 
ther to the substrate or to a previously adsorbed colloid. 
Figure [T] shows a snapshot of a typical network of con- 
nected three-patch colloids and the density profile p*(z), 
defined as the number of particles per unit volume, where 
z is the distance to the substrate. One can distinguish 
three different regimes: surface layer, liquid film, and in- 
terfacial region. As we move away from the surface, the 
density rapidly decreases (surface layer, z < 50) until a 
saturation value pi : which is constant within the liquid 
film (50 < z < 150), and vanishes in the interfacial region 
(z > 150). In the following, we discuss in detail each of 
these regimes. 

Surface layer. As the colloid-substrate interaction is 
isotropic, the patches of colloids directly adsorbed on 



FIG. 5. (color online) Analysis of the dependence on the 
fraction of two-patch colloids, rr>- (a) Liquid-film density 
(black squares) and the contribution from two- (red triangles) 
and three-patch (green spheres) colloids, (b) Total number of 
chains of two-patch colloids between three-patch ones with 
sizes 0, 1, and 2. The peak of the distribution of chains of 
unit size coincides with the peak of the film density. Results 
for 60 particle layers adsorbed on a substrate with L = 32 
averaged over 500 samples. 



the substrate are oriented randomly. By contrast, the 
colloid-colloid (i.e., patch-patch) interaction is strongly 
anisotropic and the network chains extend only along 
the direction of the patches. As the front of the network 
propagates, its branches screen the inner layers and span 
along the lateral direction, with a consequent decrease of 
the density. This effect was also reported for DLD, where 
a power-law decay with z is observed, with multifractal 
scaling p~9j[27]. However, here the conspiracy between 
the order promoted by the patches and the disorder im- 
posed by the first layer results in an exponential decrease, 
at least for the lateral sizes that were investigated. 

The plot in Fig. [2] shows the dependence of the first- 
layer density, p*(0), on the size of the substrate L and dif- 
fusion coefficient D. While there is no significant finite- 
size effect, an increase of p*(0) with D is observed (see 
snapshots in the same figure). For the irreversible adhe- 
sion of colloids on substrates, in the limit where parti- 
cles only stick to the substrate and not on top of other 
particles, extensive simulations have shown no significant 
dependence on the diffusion coefficient [28, 29 . In that 
case, the structure of the film should resemble that of 
Random Sequential Adsorption (RSA) HQ EQJ EH- In- 
stead, the patch-patch interaction promotes the forma- 
tion of colloid chains hindering the access to the sub- 
strate. Additionally, in the diffusion process, the typical 
colloid step length increases with D. Since multiple colli- 
sions can occur before irreversible binding, the larger the 
step length the higher the probability that the colloids 
can squeeze into the fjords and, eventually, arrive at the 
substrate. Consequently, the first-layer density increases 
with D towards the RSA limit (limit line in Fig.j2|. How- 
ever, patchy particle systems will always form networks 
and the RSA limit is never reached. 

Liquid film. As the lateral growth of the network pro- 
ceeds, the finite size of the substrate induces a saturation 
of the density at p = p\. For the range of values of L and 
D considered here, the network is always a fractal with 
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fractal dimension df — 2.58 ± 0.04, calculated using the 
box counting algorithm; this value is compatible with 
the one reported for DLA and DLD [19]. Notwithstand- 
ing, pi depends on L and D (see Supplementary Infor- 
mation). In both cases, a power-law scaling is observed, 
namely, p x ~ L d f~ d D~P , with d f - d = -0.4 ±0.1 and 
P — 0.23 ± 0.02. While the exponent in the size depen- 
dence results from the scale invariance of the network and 
is straightforwardly connected with df, the dynamic ex- 
ponent /3 is new and there is no equivalent in previously 
studied models. 

Interfacial region. The density vanishes at the in- 
terface between the network and the solution, which 
corresponds to the active front of the film. The po- 
sition of the interface depends on the total number 
of adsorbed particle layers N. We define Zi/ 2 as the 
thickness of the film, corresponding to the height at 
which the density is pi/ 2 (see Fig. [I]). Since the thick- 
ness of the surface layer does not depend on TV, z 1 / 2 
asymptotically scales as z 1 / 2 ~ N/ pi = NL d ~ d f . The 
profile of the density in the interfacial region scales 
as p*(z) = pi tanh [(z — £i/ 2 )6nter] , as typically observed 
for models of diffusion limited growth in the stationary 
regime (see, e.g., Ref. [33]). The decay exponent Winter 
also scales with L and D (see Supplementary Informa- 
tion). In both cases, a power law is observed, Winter ~ 
L UL D UD , with v L = -0.73 ± 0.04 and v D = 0.18 ± 0.04. 

Based on the properties of the liquid film 
and interfacial region, we propose full scaling 
of the density profile with L and D, defined as 
p*(z) = D f3 L d f~ d F [(z — z 1 / 2 )L UL D UD ] , where F [x] is 
a scaling function described by the hyperbolic tangent. 
Figure [3] shows the data collapse for several values of D 
and L. This scaling allows the definition of the network 
density for any size of the substrate and of the diffusion 
coefficient. 

Adsorption of two- and three-patch colloids. Finally we 
consider the adsorption of a mixture of two- and three- 
patch colloids. While the two-patch colloids (patches on 
the poles) favor long chains, the three-patch particles pro- 
mote branching. We define as the fraction of asorbed 
two-patch colloids and investigate how it affects the net- 
work structure. Figure [4] shows snapshots of the liquid 
film for different rjj. The larger rjj the longer the chains 
of two-patch colloids. In all cases, including vd = 0, the 
resulting network is a fractal with the same fractal di- 
mension (see Supplementary Information). However, the 
density of the liquid film can be controlled with m (see 
Fig. [5|a)). Somewhat surprisingly, a maximum is ob- 
served in the density at a fraction of two-patch particles 
around 0.35, by contrast to equilibrium coexisting gels 
(i.e. optimal networks at zero pressure) where a mono- 
tonic decrease is observed [23]. In fact, the film density 
remains above the three-patch colloid limit (rjj = 0) over 
a wide range of rp (up to around 0.6). As rp increases, 
long chains of two-patch colloids are formed and the den- 



sity is expected to decrease. However, in the absence 
of relaxation, kinetically trapped structures are obtained 
and geometrical constraints hinder the access of colloids 
to the free-patches in the inner layers. The competition 
between the formation of long chains and the maximiza- 
tion of accessible patches drives the reported film den- 
sity maximum. The plot in Fig. [5jb) shows the number 
of chains of two-patch colloids of size 0,1, and 2. With 
increasing rjj, the number of chains of size (correspond- 
ing to pairs of connected three-patch colloids) decreases 
and a maximum in the number of chains of size 1 is ob- 
served at the same concentration as the maximum in the 
film density. In the Supplementary Information we show 
that, increasing m reduces the number of free patches 
in the film and increases the average size of the chains 
exponentially. 

Conclusion. This work reveals that, in the presence 
of a substrate and under nonequilibrium conditions, new 
self-organized patterns are obtained which differ from the 
thermodynamic optimal networks or equilibrium coexist- 
ing gels (at zero pressure). For mixtures of three- and 
two-patch colloids, fractal networks of connected parti- 
cles are formed with a fractal dimension resilient over a 
wide range of diffusion coefficients and concentration of 
two-patch colloids. These networks might be of relevance 
in the fields of microfluidics and filtering. For example, 
the density variation in the surface layer resembles the 
filtering mechanism recently found in the human airways 
[34] . Here, we focused on the nonequilibrium proper- 
ties of the network, resulting from the irreversible nature 
of the binding, which dominates at low temperatures. 
However, in a significant long timescale or at higher tem- 
peratures, an adsorbed colloid might detach and rebind 
to another patch or to the substrate. In this case, the 
thermodynamic equilibrium structures might, in princi- 
ple, be reached. As a follow up, the stability and aging of 
the networks should be investigated as well as the kinetic 
pathways towards the thermodynamic equilibrium struc- 
tures. Additionally, techniques to stabilize these struc- 
tures over extended periods of time might also be a focus 
of future research. 
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SUPPLEMENTARY INFORMATION 




FIG. 6. (a) Liquid-film density for substrates with lateral size 
ranging from 16 to 128, and diffusion coefficient, D, between 
10 -3 and 1. (b) Dependence of the interfacial decay length 
Winter on L and D. 
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FIG. 7. (a) Fraction of free (unconnected) patches on two- 
patch (red triangles) and three-patch (green circles) colloids 
in the liquid film, (b) Exponential increase of the average size 
of chains of two-patch colloids. 




FIG. 8. Mass scaling as a function of the inverse box size 
in the box counting algorithm, for networks, with different 
concentrations of two-patch colloids, m- Results obtained 
with 32 layers of the 60 layers of adsorbed particles in the 
liquid film on a substrate with L = 32, averaged over 20 
samples. A fractal dimension df — 2.57 =b 0.04 was found. 



